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ABSTRACT 

We present wide area, deep, high-resolution 153 MHz GMRT observations of the NOAO Bootes field, adding to the extensive, multi- 
wavelength data of this region. The observations, data reduction, and catalogue construction and description are described here. The 
seven pointings produced a final mosaic covering 30 square degrees with a resolution of 25". The rms noise is 2 mJy beam -1 in the 
centre of the image, rising to 4-5 mJy beam -1 on the edges, with an average of 3 mJy beam -1 . Seventy-five per cent of the area has an 
rms < 4 mJy beam -1 . The extracted source catalogue contains 1289 sources detected at 5tr, of which 453 are resolved. We estimate 
the catalogue to be 92 per cent reliable and 95 per cent complete at an integrated flux density limit of 14 mJy. The flux densities and 
astrometry have been corrected for systematic errors. We calculate the differential source counts which are in good agreement with 
those in the literature and provide an important step forward in quantifying the source counts at these low frequencies and low flux 
densities. The GMRT 153 MHz sources have been matched to the 1.4 GHz NVSS and 327 MHz WENSS catalogues and spectral 
indices were derived. 
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1. Introduction 

Deep low-frequency radio surveys provide unique data which 
will help resolve many questions related to the formation 
and evolution of massive galaxies, quasars and galaxy clus- 
ters. Until now, such surveys have largely been limited by the 
corrupting influence of the ionosphere on the visibility data, 
but new techniques allow for the correction for these effects 
(e.g. |Cotton et al1|2004[ |Intema et aTT] |2009) > . Recently deep 
(0.7 - 2 mJy beam" 1 ) images have been made, in particu- 
lar with the Giant Metrewave Radio Telescope (GMRT, e.g. 
l Anant hakrishnan 2 005[) at 153 MHz ( e.g. |Ishwara-Chandra & 
Marathel 120071 ISirothia et~aT1|2009i |Ishwara-Chandra et al. 
2010). These observations can be used to study: 

Luminous radio sources at z > 4 - High redshift radio galax- 
ies (HzRGs, e.g. |Miley & De Bre uck 2008 ) provide a unique 
way to study the evolution of the most massive galaxies in the 
Universe. One of the most efficient ways of identifying these 
sources is to search for ultra-steep spectrum (USS) radio sources 
with a <, —\,S V <xV ( Rottgering et al.|1997| De Breuck et al. 
|2002| l. Low frequency observations provide an easy way of iden- 
tifying USS sources and extending these observations to lower 
flux density limits increases the distance to which these HzRGs 
can be identified. Surveying larger areas increases the probabil- 
ity of locating these rare sources. 

Distant starburst galaxies - The local radio-IR correlation 
for star forming galaxies is very tight, and seems to hold at high 
redshift (Kov acs et al.|[2006 |. However, the physical processes 
involved are poorly understood and only models that carefully 
fine-tune the time scales for the heating of the dust, the forma- 
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tion of supernovae, and the acceleration, diffusion and decay of 
the relativistic electrons can reproduce the correlation. The low- 
frequency spectral shape of galaxies reveals information about 
the amount of free-free absorption and relating this to the dust 
content, size, mass, total amount of star formation and environ- 
ment of the galaxies will further constrain the radio-IR models. 
To date, however, few galaxies have been well studied at low 
frequencies and those that have, show a diverse range of spectral 
shapes (e.g. |Clemens et al.|2Q10] >. 

Faint peaked spectrum sources - Young radio-loud AGN 
are ideal objects to study the onset and early evolution of 
classical double radio sources. They usually have synchrotron 
self-absorbed spectra and compact radio morphologies. Relative 
number statistics have indicated that these radio sources must 
be significantly more powerful at young ages, which may be 
preceded by a period of luminosity increase (e.g. Snellen et al. 
2003 ). Multi-epoch VLBI observations of individual Gigahertz 
Peaked Spectrum and Compact Symmetric Objects indicate dy- 
namical ages in the range of a few hundred to a few thousand 
years (e.g Polatidis & Conway||2003| l. Since the peak of these 
sources shifts to lower observed frequencies at higher redshift, 
low frequency observations, combined with multi-wavelength 
data, can identify these faint peakers and establish whether they 
are less luminous or at very high redshift and have different host 
properties (masses, starformation rates). 

The accretion modes of radio sources - Radio galaxies and 
radio loud quasars have been studied extensively in order to 
reveal the details of the relationship between Active Galactic 
Nuclei (AGN) and their host galaxies, in particular how their 
interaction affects their evolution. The expanding jets of radio- 
loud AGN provide a mechanism for the transfer of energy to 
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the intracluster medium and prevent the catastrophic cooling and 
formation of too-massive elliptical galaxies ( Fabian et al.|2006 



Table 1 : GMRT observation parameters for the Bootes field. 



|Best et al.|2006l[2007l[Croton et al.|2006||Bower et al.|2006| l, but 
the accretion and feedback processes and how they evolve over 
cosmic time are not fully understood. It is known that the frac- 
tion of massive galaxies which are radio-loud at z ~ 0.5 is about 
the same as observed locally (z ~ 0.1, |Best et al.||2005] l, while 
for less massive galaxies (< 10 10 5 M o ), it is an order of magni- 
tude larger. Studies of these AGN show two different types: a 
"hot" mode where radiatively inefficient accretion occurs from 
hot halo gas onto massive galaxies, and a "cold" mode where 
cold gas from major mergers drives high accretion rates. The 
strong evolution in the radio luminosity function is thus a re- 
sult of less massive galaxies experiencing more mergers and be- 
ing more active at high z. A full understanding of the differ- 
ent AGN populations, their distribution in luminosity and host 
galaxy properties, and particularly their cosmic evolution, is im- 
portant for AGN and galaxy evolutionary models. Differences in 
their host galaxy populations will provide insight into the trig- 
gering mechanisms for radio activity as well as the effect of radio 
feedback. 



In this paper we present wide, deep, high-resolution observa- 
tions of the NOAO Bootes extra-galactic field at 153 MHz taken 
with the GMRT. An initial, very deep, ~ 1 mJy beam 1 rms, 



153 MHz GMRT map of this field was presented by|Intema et al. 



( |201 l| l. Here we present additional pointings around this map ef- 
fectively tripling the size of the surveyed area at a slightly higher 
noise level. The Bootes field is part of the NOAO Deep Wide 
Field Survey (NDWFS; Jannuzi et al. 1999) and covers ~9 deg 2 
in the optical and near infra-red Bw, R, I and K bands. There is a 
wealth of additional complementary data available for this field, 
including X-ray (|Murray et al.||2005| |Kenter et aT1|2005|l, 



(GALEX;|Martin et al.|2003| l, and mid infrared (|Eisenhard t et al 



UV 



|2004l|Martin et al.|20031 l7"The region has also been surveyed at 
radio wavelengths with the WSRT at 1.4 GHz (de Vries et al. 
2002 ), the VLA at 1.4 GHz ( |Higdon et al]2005] l and 325 MHz 
( [Croft et aLl|2008| l. Recently, the AGN and Galaxy Evolution 
Survey (AGES) has provided redshifts f or 23 745 galaxies and 
AGN across 7.7 deg 2 of the Bootes field ( |Kochanek et al.|2012| l. 
This unique rich multiwavelength dataset, combined with the 
new low frequency radio data presented here, will be valuable in 
improving our understanding of the above-mentioned key topics 
in astrophysics. 



The observations presented here are the first part of the 
Two-meter Radio Mini Survey (T-RaMiSu), consisting of two 
153 MHz mosaics of similar area and depth. The second mosaic, 
centered on the galaxy cluster Abell 2256, will be presented by 
Interna et. al (in prep). 

This paper is structured as follows. In Sect.|2]we describe the 
GMRT observations of the extended region around the NOAO 
Bootes field. We describe the techniques employed to achieve 
the deepest possible images. Our data reduction relies on the 
ionospheric calibration with the SPAM package ( jlntema et al.| 
|2009| l. In Sect. [3] we describe the source detection method and 
the compilation of a source catalogue. This section also includes 
a discussion of the completeness and reliability of the catalogue 
and an analysis of the quality of the catalogue. The source counts 
and spectral index distributions are presented in Sect. [4] Finally, 
Sect.|5]summarises and concludes this work. 



Parameter 


Central 


Flanking 


Observation Dates 


3 June 200: 


3-6 June 2006 


Pointings 


Bootes 


Bootes A-F 


Primary Calibrator 


3C48 


3C48 


Total Time on Calibrator 


20 min 


51 min 


Secondary Calibrator 


3C286 


3C 286 


Total Time on Calibrator 


9 x 10 min 


10 x 4.5 min (per day) 


Cadence 


50 min 


30 min 


Total Time on Target 


359 min 


205 min (per pointing) 


Integration time 




16.9 s 


Polarisations 




RR,LL 


Channels 




128 


Channel Width 




62.5 kHz 


Total Bandwidth 




8.0 MHz 


Central Frequency 




153 MHz 


Table 2: Pointing centres of the Bootes central and flanking 


fields. 






Field 


RA 


DEC 




(J2000) 


(J2000) 


Bootes 


14:32:05.75 


+34:16:47.5 


Bootes A 


14:32:05.75 


+36:06:47.5 


Bootes B 


14:24:19.53 


+35:10:52.5 


Bootes C 


14:24:29.58 


+33:20:54.5 


Bootes D 


14:32:05.75 


+32:26:47.5 


Bootes E 


14:39:41.92 


+33:20:54.5 


Bootes F 


14:39:51.97 


+35:10:52.5 



2. Observations and Data Reduction 

2.1. Observations 

The central Bootes field was previously observed with the 
GMRT from 3-4 June 2005 ( |Intema et al-pHl) . We use the 
data from a single day of this observing run, combined with new 
observations of six flanking fields taken during 3-6 June 2006 
with the GMRT at 153 MHz. Data from the first day only of the 
first observing run, 3 June, was used as the RFI situation was 
marginally better on this day and the length of a single day's 
observation, 359 min, compares well with that of the new ob- 
servations of the flanking fields, 205 min, which leads to a more 
uniform mosaic. Table[T]lists the observational parameters used, 
highlighting any differences between the two sets of observa- 
tions. The flanking fields are arranged on a hexagonal grid with 
a radius of 110' just beyond the half power point of the primary 
beam of the GMRT at 153 MHz (6» FWHM ~ 3°); Table|2]gives the 
central coordinates of each pointing. Typically 26 - 27 of the 30 
antennas were available during each observing run. 3C48 and 
3C 286 were observed as phase, bandpass and flux density cali- 
brators. For each of the four days, the target fields were observed 
in sets of ~ 4.5 min each, followed by a calibrator observation 
(3C 286) of ~ 4.5 min. 3C 48 was observed at the beginning or 
end of each day for ~ 20 - 30 min. The frequent (~ 30 min 
interval) calibrator observations of 3C286 provide a means to 
track changes in the GMRT system, RFI and ionospheric condi- 
tions, and flux density scale. The short target field observations 
spread over each night of observing provides fairly uniform uv- 
coverage. 
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2.2. Data Reduction 

The data for the central pointing was re-reduced in the same 
manner as the new flanking fields in order to allow for consistent 
integration into a single mosaic. The data reduction consisted 
of two stages: "traditional calibration" followed by directional- 
dependent ionospheric phase calibration, both of which were 



Table 3: Final CLEANing (top) and SPAM (bottom) parameters 
for individual Bootes fields. 



implemented in Python using the ParselTongue ( Kettenis et al. 
2006) interface to the Astronomical Image Processing System 
package (AIPS; Greisen] [1998[ |. Ionospheric calibration was 
done with the "Source Peeling and Atmospheric Modelling" 
ParselTongue-based Python module (SPAM; Interna et al. 2009). 

The data for each day were calibrated separately. The flux 
density scale was set and initial amplitude, phase and bandpass 
calibration were done using 3C 48. 3C 48 is brighter than 3C 286 
and provides a better determination of the bandpass. In order to 
reduce the data volume, the LL and RR polarisations were com- 
bined as Stokes I and every 5 channels were combined to form 
18 channels of width 0.3125 MHz yielding an effective band- 
width of 5.625 MHz. After this calibration, the Mv-data from all 
four days for each target, were combined. 

Initial imaging of each target field was done after a phase- 
only calibration against a model field constructed from NVSS 
sources within each field. Table [3] lists the important imag- 
ing parameters. The calibration was then improved by several 
rounds of phase-only self-calibration followed by one round of 
amplitude and phase self-calibration where gain solutions were 
determined on a longer time-scale than the phase-only solu- 
tions. Excessive visibilities were determined from the model- 
subtracted data and were removed. Additional automated re- 
moval of bad data causing ripples in the image plane was done 
by Fourier transforming the model-subtracted images and iden- 
tifying and removing extraneous peaks in the Mv-plane. Further, 
persistent RFI was flagged and low level RFI mo delled and sub - 
tracted using the LowFRFif] routine in ObitTalk ( |Cotton|2 008). 
After self-calibration the rms noise in the inner half of the pri- 
mary beam area was 2.5 mJy beam -1 in the central field and 
3.5 - 5 mJy beam" 1 in the flanking fields, with the local noise 
increasing 2-3 times near the brightest sources. Note, the pres- 
ence of extremely bright sources with peak flux densities of the 
order of 5 - 8 Jy beam" 1 prior to primary beam correction in 
flanking fields D through F resulted in the slightly higher overall 
noise in these fields. 

Significant artefacts, however, remained in all fields near 
bright sources. To reduce these we applied the SPAM algorithm 
on the self-calibrated data. The SPAM parameters are listed in 
the bottom part of Table [3] which include the number of iono- 
spheric layers modeled and their heights and relative weights, 
the slope of the assumed power law function of phase structure 
resulting from turbulence (y) and the number of free parameters 
in the fit ; seellntema et al. (2009) for a more detailed description 
of the meaning of these parameters. Three iterations of peeling 
were done: in the first we only applied the peeling solutions to 
the peeled sources and in the final two we fitted an ionospheric 
phase screen to the peeling solutions. Up to 20 sources with flux 
densities above 0.4 Jy (not corrected for primary beam effects) 
were peeled in the final stage in each field. The screen was made 
up of two equally-weighted turbulent layers at 250 and 350 km. 
SPAM also allowed for the determination of and correction for 
antenna-based phase discontinuities. 

In order to have a homogeneous point spread functions in all 
pointings, final images were made with a circular restoring beam 



Parameter 


Value 


Widefield imaging 


polyhedron facet-based" 




multi-frequency synthesis 6 


Deconvolution 


Cotton-Schwab Clean c 


field size 


4° 


facets 


85 


facet size 


32.4' 


facet separation 


26.4' 


Weighting 


Robust' -0.5 d 


UVBXFN, UVBOX 


4, 1 


Clean box threshold 


5 cr 


Clean depth 


3o- 


Pixel size 


3.8" 


Restoring beam 


25" circular^ 


Spam calibration cycles 


3 


Peeled sources 


20 7 


Layer heights (weights) 


250 km (0.5) 




350 km (0.5) 


Turbulence parameter y 


5/3* 



Model parameters 
Reference catalogue 



< 20 
NVSS'' 



Perley[) 1989); Cornwell & Perley|l |1992j 



Conway et al. ( 1990) 



(c) [Schwab| |l9P| l^ TCoTton| ( | 1 999 1 



Final imaging parameters 
(e) |Briggs|jl995) 
(r) 14 for field F 
te) pure Kolmogorov turbulen ce 
w |Condon et al.N 199411 1998 1 



Cornwell et al. 



1999) 



of radius 25" and a pixel size of 3.8". The flux density scales of 
the central and flanking fields were scaled up by 60 per cent 
and 30 per cent respectively based on information from 3C 286 
(discussed in Sect. |3.4) . In the final individual field images, the 
rms noise in the central half of the primary beam area before 
primary beam correction was 1.8 mJy beam" 1 in the central field 
and 2.5-2.7 mJy beam" 1 respectively in the flanking fields. This 
is 3 - 5 times the theoretical noise, similar to the factor above 



1 Obit Development Memo Series # 16 see http : //www . cv . nrao . 
edu/-bcotton/Obit.html 



thermal noise obtained by deeper single pointing of Interna et al. 
( |201 1) . The seven pointings were each corrected for the primary 
beam of the GMRT up to a radius of 1.6°, where the primary 
beam correction factor drops to 40 per cent of its central value, 
and were then mosaicked together by weighting the final image 
by the inverse of the square of the rms noise of each individual 
pointing. Figure[T]illustrates the variation in rms noise across the 
mosaic which is shown in entirety in Fig. [2] The noise level is 
smooth and around 2 mJy beam" 1 across the interior of the map, 
and increases towards the edges to about 4-5 mJy beam" 1 . The 
average noise in the final mosaic is 3.0 mJy beam" 1 , with 49 per 
cent under 3 mJy beam" 1 and 74 per cent under 4 mJy beam" 1 . 
A small portion of the mosaic covering the inner square degree is 
shown in Fig.[3]to illustrate the resolution and quality of the map. 
There remain some phase artefacts visible around the brightest 
sources, which have not been entirely removed during peeling. It 
is possible that some artefacts are caused by elevation-dependent 
pointing errors, since each pointing was observed in a series of 
scans with varying elevations (Tasse et al. 2007 Mohan et al. 
[200T] [Chandra et al.|2004| >. 
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21S.00 
RA 1)2000) 



Fig. 2: Greyscale map showing the entire mosaic. The image covers 30 square degrees. The greyscale shows the flux density from 
-3(T avg to 25cr m ,g where cr avg = 3.0 mJy beam -1 is the average rms across the entire mosaic. 



3. Source Detection and Characterisation 

3.1. Detection 

We used the PyBDSM packag^] to detect and characterise 
sources in the mosaic image. PyBDSM identifies islands of con- 
tiguous emission by identifying all pixels greater than the pixel 
threshold and adding each of these pixels to an island of con- 
tiguous pixels exceeding the island threshold. Each island is fit 
with one or more Gaussians which are subsequently grouped 
into sources. Sources are classified as 'S' for single sources, 
'M' for multiple-Gaussian sources and 'C' for components of a 
multi-source island. From the fitted parameters the deconvolved 
sizes are computed assuming the theoretical beam. Errors on the 
fitted parameters are computed following |Condon| ( | 1 997} . Prior 
to source detection the local background rms is determined by 



measuring the pixel statistics within a sliding box. For determin- 
ing the rms background in our map we used a box size of 100 
pixels to capture the variation in local noise around the brightest 
sources. We used a pixel threshold of 5cr L and an island thresh- 
old of 3cr L . In generating a source list we allowed all Gaussians 
in each island to be grouped into a single source. PyBDSM de- 
tected 1296 sources from 1578 Gaussians fitted to 1301 islands, 
of which 1073 were single-component 'S' sources. Based on vi- 
sual inspection a small number of sources were removed as they 
were false, or bad, detections on the edge of the image. 

The final catalogue consists of 1289 sources between 4.1 
mJy and 7.3 Jy and is available as part of the online version 
of this article and from the CDS0 The flux scales of the indi- 
vidual pointings were adjusted prior to mosaicing as described 
in Sect. 3.4 and the astrometry in the catalogue has been cor- 



http : / /home . strw . leidenuniv . nl/~mohan/anaamika 



http : // cdsweb . u- strasborg . fr/ 
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0.04 g 



218.00" 
RA 1)2000) 



Fig. 3: Zoom-in of the central part of the mosaic. The image covers 1 square degree. The greyscale shows the flux density from 
-3cr m , g to 25o- avg where cr avg = 3.0 mJy beam -1 is the average rms across the entire mosaic. 



rected for a systematic offset also described in Sect. 3.4 A sam- 
ple of the catalogue is shown in Table [4] where the columns are: 
(1) Source name, (2,3) flux- weighted position right ascension, 
RA, and uncertainty, (4,5) flux-weighted position declination, 
DEC, and uncertainty, (6) integrated source flux density and un- 
certainty, (7) peak flux density and uncertainty, (8,9,10) fitted 
parameters: deconvolved major- and minor-axes, and position 
angle, for extended sources, (11) local rms noise, and (12) the 
number of Gaussians fitted to the source. Extended sources are 
classified as such based on the ratio between the integrated and 
peak flux densities (see Sect. 3.2 1. Unresolved sources have a 
'-' listed for all their fitted shape parameters (semi-major and - 
minor axes and position angle) or for only the semi-minor axis 
where the source is resolved in one direction. For extended 
sources consisting of multiple Gaussians, the fitted parameters 
for each Gaussian are given on separate lines in the table, listed 



as 'a', 'b', etc. Images of the 25 brightest sources are shown in 
Appendix [A] 

3.2. Resolved Sources 

In the presence of no noise, the extendedness of a source can 
simply be determined from the ratio of the integrated flux density 
to the peak flux density, S i/S p > 1. However, since the errors on 
Si and S p are correlated, the S ,/S p distribution is skewed, par- 
ticularly at low signal-to-noise. To determine an upper envelope 
of this distribution, we performed a Monte-Carlo simulation in 
which we generated 25 random fields containing ~ 10000 ran- 
domly positioned point sources with peak flux densities between 
0.1 cr and 20cr, where cr was taken to be 3 mJy beam -1 . The 
source flux densities are drawn ra ndom ly from the source count 



distribution, dN/dS oc S (Sect. 4.1 1. We neglect the deviation 



of the true source counts from a power law slope at high fluxes 
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Table 4: Sample of the GMRT 153 MHz source and Gaussian-component catalogue. 



Source ID 




0" RA 


1 J 1 A 


0" £>£C 


C . 
^ / 


c 

b l' 


Li 


D 


<P 


Tffl s 


AT b 

''gauss 




[deg] 


r"i 


[deg] 


["] 


[mjy] 


[mjy beam~^] 


r" l 
1 J 


1 1 


[deg] 


[mjy beam - 1 ] 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Tl A AT} t i "J C f\H 

J 144733+350/ 


221.85/96 


0.9 


35.13126 


1.5 


158 ± 33 


69 ± 15 


38.3 ± 3.4 


19.3 ±1.9 


5 ± 6 


5.4 


1 


T 1 A A £lCO i nAO 

J 144658+3308 


221.74358 


2.9 


33.14058 


1 . 1 


42 ± 13 


20 ± 6 








4.9 


2 


a 


111 "7/10/1") 

221.74843 


2.8 


33.14096 


1 .8 


19 ± 8 


19 ± 5 












b 


221.73914 


2.8 


1 1 A f\f\ 1 

33.14001 


2.6 


23 ± 8 


18 ± 5 


15.9 ± 6.8 


9.5 ± 5.6 


56 ± 90 






T1 A A*lf\C i 1 A A^ 

J 144705+3442 


221.77266 


1.8 


34.71605 


1 .6 


32 ±9 


21 ±1 








4.7 


1 


t 1 A A"~i r\/i , Tin 

J 144706+3457 


221.77552 


1.4 


34.95096 


1.5 


25 ± 9 


28 ± 7 








5.0 


1 


T l A A £L A C i 

J 144645+3330 


221.69153 


1.8 


33.50822 


1 .4 


31 ± 9 


29 ± 7 








4.8 


1 


Tl A A £ A f\ i OO^^ 

J 144640+3322 


221.66919 


0.9 


33.36920 


0.8 


72 ± 16 


62 ± 13 








4.9 


1 


T1 A A £. A £. i A Af\ 

J 144646+3440 


221.69543 


1.5 


34.67482 


1 . 1 


29 ±9 


32 ± 8 








5. 1 


1 


J 144648+3546 


221. 70230 


1.9 


t c n /i 1 1 

35. 77421 


1 .4 


19 ± 8 


25 ± 6 








5.0 


1 


t i A a 1 ^ i ")ir\"> 

J 144613+3303 


111 CCOT 1 

221.55831 


0.7 


33.06494 


0.6 


169 ± 35 


131 ± 27 








5.0 




J 144638+3553 


221.65996 


2. 1 


35.88845 


2.3 


37 ±11 


27 ± 7 








5.7 


1 


J 144626+35 12 


221.61 124 


1.0 


35.2092 7 


0.7 


314 ± 65 


244 ± 50 








7.5 




J 1 '\'\\J 1 7TjtZJ 


99 l SRI 9^ 


3 9 


^4. A1AA1 
JM-.M-ZH-H / 


2 4 


46+12 


26 + 7 










y 


J144606+3316 


221.52865 


2.9 


33.26922 


2.2 


20 ±7 


17±5 








4.2 


1 


J 144557+3251 


221.49062 


0.6 


32.86032 


0.5 


122 ±26 


111 ±23 








5.0 


1 


J 144555+3237 


221.48098 


2.4 


32.62328 


3.6 


24 ±7 


17 + 5 








4.4 


1 


J144617+3506 


221.57405 


1.1 


35.11639 


1.3 


158 ±34 


75 ± 16 








5.3 


2 


a 


221.57091 


0.9 


35.11314 


0.7 


93 ±20 


78 ± 16 


15.3 ± 1.8 


5.6 ± 1.4 


72 ± 16 






b 


221.57923 


1.4 


35.12223 


2.1 


65 ± 15 


38 ±9 


29.4±5.1 


10.6 ±2.8 


12± 14 






J144602+3339 


221.50949 


2.1 


33.66347 


1.3 


42 ± 10 


27±6 


29.3 ±5.3 




64 ±11 


3.8 


1 


J144607+3503 


221.52967 


1.2 


35.05975 


0.6 


81 ± 18 


62 ± 13 






178 ±6 


4.8 


1 



Parameters are given for extended sources to which Gaussian components were successfully fit. 

A "-" indicates a poor Gaussian fit. In these cases the total flux density quoted is the total flux density in the source island. 




RA (J2000} 

Fig. 1: Greyscale map showing the local rms noise measured 
in the mosaic image. The greyscale shows the rms noise from 
0.50-^ to 2<r avg , where cr m , g = 3.0 mjy beam -1 is the ap- 
proximate rms in the mosaic centre. The contours are plotted 
at [1/V2.1, V2] X (Tavg- Peaks in the local noise coincide with 
the locations of bright sources. 



as there are very few sources at these fluxes. The rms noise map 
for these fields was taken from the central 4000 x 4000 pixel 2 
of the residual mosaic. Source detection was performed in the 
same manner described in Sect. |3.1| thus only ~ 750 sources 
in each field satisfy the detection criterion of peak flux density 
> 5cr. The S j/S p distribution produced from the Monte-Carlo 
simulation is plotted in the left panel of Fig. [4] To determine the 



95 per cent envelope, a curve was fit to the 95th percentile of 20 
logarithmic bins across signal-to-noise ratio. The fitted envelope 
is characterised by: 

Si/S p = 1 + |(0.01 + 0.02) 2 + (3.58 ± 0.10) 2 (cr L /S pf^ . 

The measured distribution of S , fS p as a function of signal- 
to-noise ratio is shown in the right panel of Fig. |4] The line 
shows the upper envelope from the Monte-Carlo simulation. Of 
the 453 sources that lie above this line (35 per cent of all 1289 
sources), approximately 41, i.e. 9 per cent, are not truly extended 
and merely lie above the line by chance. However, all these 
sources above the line are listed in the catalogue as extended 
and the measured deconvolved FWHM major and minor axes 
are given. 

3.3. Completeness and Reliability 

To quantify the completeness and reliability of the catalogue, we 
performed a similar Monte-Carlo simulation to that described in 
the previous section. However, in this case approximately 25 per 
cent of the artificial sources inserted into the noise map were ex- 
tended sources - Gaussians with FWHM larger than the beam- 
size. This allows for a better estimate of the completeness and 
reliability in terms of integrated flux densities. 

The completeness of a catalogue represents the probability 
that all sources above a given flux density are detected. We have 
estimated this by plotting the fraction of detected sources in our 
MC simulation as a function of integrated flux density (left panel 
of Fig. [5j, i.e. the fraction of input sources that have a cata- 
logued flux density using the same detection parameters. Due 
to the variation in rms across the image, the detection fraction 
has been multiplied by the fraction of the total 30 deg 2 area in 
which the source can be detected. The completeness at a given 
flux density is determined by integrating the detected fraction 
upwards from a given flux density limit and is plotted as a func- 
tion of integrated flux density in the right panel of Fig. [5] We 
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Fig. 4: Left The simulated ratio of integrated to peak flux density as a function of signal-to-noise ratio for sources from the 25 Monte- 
Carlo simulations. For 20 logarithmic bins in signal-to-noise ratio, the black points show the threshold below which 95 per cent 
of the sources lie in that bin. The red line shows a fit to this upper envelope. Right The measured ratio of integrated to peak flux 
density as a function of signal-to-noise ratio. The line shows the upper envelope containing 95 per cent of the unresolved sources as 
determined from Monte -Carlo simulations. 




0.0 1 ■ ■ — — ' ■ ■ ■ ■ ■ ■ — — 1 0.0 1 ■ ' — — ' ' ■ ■ ' ' ' — — 1 

I0 1 10 2 I0 1 10 2 

S, [mjy] S, [mjy] 

Fig. 5: Left Fraction of sources detected as a function of integrated flux density to local noise ratio calculated from 25 Monte-Carlo 
simulations. The solid line shows the mean of all 25 randomly generated fields and the two dotted lines show the lcr uncertainty. 
Right Estimated completeness of the catalogue as a function of integrated flux density limit accounting for the varying sensitivity 
across the field of view. 



thus estimate that the catalogue is 95 per cent complete above a 
peak flux density of 14 mjy. 



The reliability of the catalogue indicates the probability that 
all sources above a given flux density are real. In the left panel 
of Fig. [6] the false detection rate FDR, i.e. the fraction of cat- 
alogued sources that do not have an input source, is plotted as 
a function of the integrated flux density. Integrating up from a 
given detection limit and multiplying by the normalised source 
flux distribution, we can determine an estimate of the overall 
FDR or reliability, R = 1 — FDR, of the catalogue. The relia- 
bility is plotted as a function of integrated flux density limit in 
the right panel of Fig. [6] For a 14 mjy detection threshold, the 
reliability is 92 per cent. 



3.4. Astrometric and Flux Uncertainties 

Errors in the phase calibration introduce uncertainties in the 
source positions. To assess these uncertainties and determine any 
systematic offsets we selected a sample of sources with peak 
flux densities at leas t 10cr £ . We searched for 1.4 GHz NVSS 
dCondon et al.||1998b sources within 45" of these targets. 745 



matches were found. From this sample, we measured a small 
offset of (Aa, A6) = (0.44", -0.21"), which is of the order of the 
pixel size of the 153 MHz observations and the NVSS accuracy 
(~ 1"). A correction for this offset has been applied to all sources 
in the catalogue. The scatter in the offsets between the GMRT 
and NVSS positions is a combination of noise-independent cal- 
ibration errors, e, in both the GMRT and NVSS data as well 
as a noise-dependent error, <x, from position determination via 
Gaussian-fitting: 



(r = e 



GMRT 



NVSS 



GMRT 



+ cr- 



NVSS 



1 



W. L. Williams et al.: T-RaMiSu: The Two-meter Radio Mini Survey. I. 




S, [mjy] S, [mjy] 

Fig. 6: Left False detection rate as a function of peak flux density to local signal-to-noise ratio calculated from 25 Monte-Carlo 
simulations. The solid line shows the mean of all 25 randomly generated fields and the two dotted lines show the lcr uncertainty. 
Right Estimated reliability of the catalogue as a function of integrated flux density limit accounting for the varying sensitivity across 
the field of view. 



From |Condon et al. ( |1998| >, the NVSS calibration errors are 
(fa, £s)nvss — (0.45", 0.56")- To separate the noise-dependent 
and -independent uncertainties we select from the above sam- 
ple only the NVSS sources with position errors of less than 0.6" 
and measure a scatter of (cr , <t$)cmrt = (0.67", 0.65"). For this 
very high signal-to-noise sub-sample of 107 sources the noise- 
dependent fit errors for both the GMRT and NVSS can safely 
be assumed to be small so we determine the GMRT calibra- 
tion errors to be (e a ,e$)GMRT = (0.50", 0.32"). These are added 
quadratically to the Gaussian-fit position uncertainties in the cat- 
alogue. 

Similarly, in addition to the noise-dependent Gaussian fitting 
uncertainties on the fluxes ( |Condon 1997| l, the uncertainty in the 
measured flux densities all consists of a noise-independent com- 
ponent. The uncertainty introduced through transferring the flux 
density scale from the calibrator to the target fields is the main 
such uncertainty and depends on a number of factors: (i) the data 
quality, (ii) the accuracy of the model, and (iii) differences in ob- 
serving conditions between the calibrator and target. 

Like the target data, the calibrator data is adversely affected 
by RFI and the ionosphere. The RFI conditions of the flanking 
field observations were similar to those when the central point- 
ing data were taken, however, the ionosphere was not as calm. 
Following Interna et al. ( 201 1\ we adopt a slightly inflated, ad- 
hoc amplitude uncertainty of ~ 4 per cent due to RFI and iono- 
spheric effects. 

The calibrator model is of a point source whose flux density 
at 153 MHz is predicted from the Perley-Taylor model based 
on flux density measurements at many frequencies. 3C48 is a 
point source of 64.4 Jy at 153 MHz. [lntema et al.] ( |201 1} provide 
an improved model for 3C286, a point source of 31.01 Jy at 
153 MHz, and estimate a flux density uncertainty of 5 per cent. 
The large field of view, however, means that there are other 
fainter sources present in the calibrator field. For similar dura- 
tion observations of 3C 286 [Tntema et aL (2011 1 set an upper 



limit of 1 per cent on the flux density uncertainty. Since 3C 48 is 
about a factor of two brighter, we estimate that the flux density 
uncertainty due to additional sources in the 3C 48 field is also at 
most 1 per cent. 

Individual antennas are sensitive to the galactic diffuse radio 
emission which varies across the sky and so may be different 



for the calibrator and target fields thereby introducing an off- 
set to the flux density scale as well as additional uncertainty. 
However, since the GMRT does not measure the sky tempera- 
ture, we require external information to take this into account. 
Following |Tasse et al. (|2007) and Interna et al. ( 201 1\ we de- 
termine the mean off-source sky-temperature from the Haslam 
al.| ( |1982| l all-sky radio maps at 408 MHz: both the Bootes 
and 3C 286 fields have sky temperatures of ~ 20 + 1 K. Using 
the equation from |Tasse et al.| ( [2007| l, this implies that no offset 
in the flux density scale is required for 3C 286 and we estimate 
a gain uncertainty of 2 per cent. However, the sky temperature 
near the primary calibrator 3C 48 is 24+ 1 K which implies a flux 
density correction of 0.92 with an estimated uncertainty of 8 per 
cent. Since the flux density scale is linear, this offset is applied 
post hoc to the measured flux densities. 



Prior to combining the individual pointings, we compared 
the measured primary beam-corrected flux densities of sources 
in the overlapping regions (approximately 110 - 150 sources 
per region) and found those in the flanking fields to be consis- 
tently higher by 30 + 5 per cent. To investigate this we made 
images after calibration using 3C 286 as the primary calibrator. 
This yielded consistent fluxes between the central and flanking 
fields. It is likely that significant time-dependent changes in the 
GMRT systems over the course of each observing night were 
captured by the regular (each 30 min) observations of 3C 286. 
We thus used the 3C286-calibrated images to derive a correc- 
tion to the flux density scales of the flanking fields, a factor of 
1.3, before combining the individual pointings. The uncertainty 
of this correction is 10 per cent. 



The total estimated uncertainty in transferring the flux den- 
sity scale is of the order of 20 per cent which we add quadrat- 
ically to the measured Gaussian fit uncertainty for each source. 
Comparison of the flux density of bright sources measured in 
the individual pointings after the above correction shows good 
agreement between the flux density scales of the individual 
pointings and the measured scatter is ~ 16 per cent, which also 
includes a contribution by the noise-dependent terms. 
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-1 1 1 1 1 r 

GMRT 01 
J142132 + 3513 
2.985 mjy/beam 
RA = 215.385 
DEC = 35.220' 



of the area in which it can be detected (e.g. Wind horst et al.| 



215.52' 215.40° 215.28° 

RA (J2000) 



> 4 

1- 



i GMRT D2 
X[lj41pS + 3401 
r = 2.621 it#/beam 
RA = 220.485 
DEC = 34.026 



j2l 



220.68° 220.56° 220.44° 

RA (J2000) 

Fig. 7: Postage stamps showing Dl, RA=14:21:32, 
DEC=+35:12:12(fop) and D2, RA=14:41:56, DEC=+34:01:34 
(bottom). The greyscale goes from 0.5cr to 5<x and the images 
have been smoothed with a Gaussian of 50 arcsec. WENSS 
contours are overlaid at [1.5,3.0,10.0] x <x L where cr L is 
the local rms in the WENSS images - 3.5 mJy beam 1 and 
3.7 mJy beam 1 respectively for Dl and D2. 



3.5. Diffuse Sources 

We have identified two faint diffuse sources in the final mosaic 
which were not detected by PyBDSM as their peak flux densi- 
ties are too low. Postage stamps of these two sources are shown 
in Fig. [7] The first, Dl, is located at RA = 14:21:32, DEC = 
+35:12:12. This source has previously detected in WENSS by 
Delain & Rudnick (2006 ) who have associated it with a galaxy 
group at z — 0.01. The second diffuse source, D2, is located at 
RA = 14:41:56, DEC = +34:01:34. 



4. Analysis 

The 1289 sources in the catalogue provide a statistically signif- 
icant sample across three orders of magnitude in flux density 
from 4 mJy to 7 Jy. In this section we present the derived 153 
MHz source counts and spectral index distributions based on 
matching these sources to catalogues at 1.4 GHz and 320 MHz. 



4.1. Source Counts 

The Euclidean-normalized differential source counts are shown 
in Fig. [8] Due to the large variation in rms across the mosaic, the 
sources are not uniformly detected across the image, i.e. faint 
sources can only be detected in a smaller area in the inner part 
of the image. We therefore weight each source by the inverse 



1985 ), which also accounts for the varying detection area within 



a given flux density bin. Accurate derivation of the source counts 
is complicated by a number of effects. In general, noise can scat- 
ter sources into adjacent bins, most noticably at low flux densi- 
ties. A positive bias is introduced by the enhancement of weak 
sources by random noise peaks (Eddington bias). Furthermore, 
low surface brightness extended sources can be missed as their 
peak flux densities fall below the detection limit. We have used 
our Monte-Carlo simulations to estimate the combined contribu- 
tion of these effects and derive a correction factor to the observed 
source counts. Errors on the final normalised source counts are 
propagated from the errors on the correction factors and the 
Poisson errors (Gehrels 1986 1 on the raw counts per bin. The 
flux density bins start at three times the average rms, 15 mJy, 
and increase in factors of 2^ 4 , 2^ 2 or 2 chosen to provide source 
counts of 60 - 80 in most, except for the highest, flux density 
bins. Table [5] lists (i) the flux density bins, (ii) the central flux 
density of the bin, (iii) the raw counts, (iv) the effective detection 
areas for sources at the lower and upper limits of the flux density 
bin, (v) the effective area corresponding to the bin centre, (vi) 
the mean weight of the sources in the bin, (vii) the correction 
factor, and (viii) the corrected normalised source counts. 

We have compared our source counts with the little observa- 
tional data available at this frequency. Our source counts agree 
well with those derived by Interna et al. ( 201 l| l for the cen- 
tral field. Since their image is approximately three times deeper 
than our mosaic, the good agreement at low flux densities lends 
credance to our correction factors. The recent source counts from 

20T0) l 



Ghosh et al. ( 2012)) and those by llshwara-Chandra et al. 



for a smaller, slightly shallower GMRT field also agree well with 
our data, except the Gho sh et aI7| ( |2012} ) counts deviate at low flux 
densities, becoming increasingly lower. At the high flux end, the 
7C 151 MHz source counts (McGilch rist et al.|1990] > match our 
counts well. We have fit a power law over the flux density range 
15-400 mJy and obtain a slope of 0.93+0.04 which is consistent 



with, but slightly steeper than, the 0.91 obtained by Interna et al. 
( |201 1) across the same flux density range. Likewise, it is con- 



sistent with the value of 1.01 found by Ishwara-Chandra et al. 



(2010), but is slightly shallower. The source counts derived from 
the small sample o f|George & Stev ens (2008 ) (not plotted) are fit 
by a single power law with a slope of 0.72, but their deviation is 
probably due to poor statistics. Model source counts have been 
derived by |Wilman et al.| ( [2008| ) for the 151 MHz source pop- 
ulation predicted from the extrapolated radio luminosity func- 
tions of different radio sources in a ACDM framework. The 
[Wilman et al.| ( |2008| ) model catalogue has been corrected with 
their recommended post-processing, which effectively reduces 
the source count slightly at low flux densities. The dominant 
source population at flux densities above ~ 200 mJy is that of 
FRII radio sources. Only below this flux density does the FRI 
population begin to dominate. There is a general agreement be- 
tween our data and this model which has an approximate power- 
law slope of 0.79 between 10 and 400 mJy. At low flux densities 
it is likely that the Wilman et al. (2008 ) counts slightly overesti- 
mate the true counts due to double counting of hybrid AGN-star 
forming galaxies. 

4.2. Spectral Index Distributions 



While deep 1 .4 GHz data exists for the Bootes Field ( |de Vries] 
et al.||2002 i, this only covers the central 7 deg 2 . This data was 



used in Interna et al. ( 201 1\ in a 153 MHz flux-limited spec 
tral index analysis. However, we choose to compare our source 
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Fig. 8: Euclidean-normalized differential source counts for the GMRT 153 MHz catalogue (filled black circles) in 18 logarithmic 



flux density bins between 15 mJy and 6.5 Jy. For comparison we have plotted the 153 MHz source counts from Interna et al. 
(open squares) for the central Bootes pointing, from Ghosh et al. (|2012|l (blue filled squares) and from Ishwara-Chandra et al 



(2011 



(2010: 



(red filled triangles), as well as the 151 MHz source counts from McGilchrist et al. ( 1990| > for part of the 7C catalogue (open inverted 
triangles). Also shown is a source count model by Wilma net al.| ( |2008| ) (dotted line with shaded area indicating the lcr errors) and 
a power law fitted between 150 - 400 mJy (dashed line) which has a slope of 0.93 + 0.04. 



list to the NVSS 1.4 GHz catalogue (Cond on etal.|1994| > which 
covers our entire survey area at a comparable resolution. We 
searched for NVSS counterparts within 45" of each GMRT 
source. Despite the relatively small difference in resolution be- 
tween the NVSS (45") and the GMRT (25") data, a small num- 
ber of GMRT sources (9 pairs) were matched the same NVSS 
source. Also, due to differences in the grouping of compo- 
nents into sources, we merged 16 pairs of NVSS sources which 
matched a single GMRT source. Sources were merged by sum- 
ming their total flux densities. A spectral index was calculated 
for each GMRT source based on the combined flux density of 
merged sources. 

We matched 1 134 NVSS sources to 1 127 GMRT sources and 
then used this matched subsample to compute the spectral inde?^ 
distribution which is shown in Fig. [9] The flux density limit of 
2.5 mJy at 1.4 GHz biases the detection of 1.4 GHz counter- 
parts to fainter 153 MHz sources to those with flatter spectra. 
168 GMRT sources have no match in NVSS. These are consis- 
tent with having steeper spectral indices below the diagonal line 
in Fig. 9]and we therefore provide an an upper limit to the spec- 
tral index given the NVSS flux density limit. The mean spectral 
index is -0.87 + 0.01, calculated using the Kaplan-Meier esti- 
mator (KM; e.g. Feigels on & Nelson||1985 1 to account for the 



upper limits. This value is comparable to those found by |Intema 
et^([20rTl), -0.79, [Ishwara-Chandra & Marathe|(|2007l), -0.85, 



Sirothia et aE1 < |2009| >, -0.82, and Ishwara-Chandra et aT](|20T0}, 
-0.78. By considering the KM mean spectral index within 5 



logarithmic flux density bins between 85 mJy and 1 Jy (over- 
plotted in Fig. [9] and listed in Table [6]), we find a gradual steep- 
ening of the spectral index with increasing flux density, from 
~ -0.84 at ~ 30 mJy to ~ -0.97 at F > 600 mJy. This trend 
is still clear if the first flux density bin is ignored (i.e. consider- 
ing F > 40 mJy) assuming that this bin remains biased by the 
upper limits. This is consistent with what is found in the litera- 
ture (e.g. Ishwara-Chandra et al.|2010) Tasse et al.|2006 Cohen 



|et al.|2004||de Vries et al. 2002 ). Since it appears that there is no 
spectral steepening or flattening due to redshifted curved spectra 
(Bornancini et al. 2010), this flattening is likely due to a cor- 
relation between source luminosity and spectral index (P - a), 
which is known to exist for FRII radio ga laxies (e.g. |Blundell 



let al. 



1999). According to the models of Wilma n~eTal.| ( |20081 >, 



The spectral index is denned as 5 „ oc V 



the observed 153 MHz source population is dominated by FRII 
galaxies at these flux density levels (> 20 mJy). 

We also compared our source list to WENSS at 327 MHz 
( Rengelink et al.|1997| l, noting that the errors in this spectral in- 
dex are much greater due to the smaller difference in frequency. 
The WENSS beam is 54" x 54"/ sin 6, or 54" x 96" at the dec- 
lination of the Bootes field. We thus searched for WENSS coun- 
terparts within 96" of each GMRT source. Of the 1289 GMRT 
sources we matched 689 to 675 WENSS sources. The 14 pairs of 
GMRT sources within the beam of a single WENSS source were 
combined as described in the previous paragraph and spectral 
indices determined for each based on the combined flux den- 
sity. A visual check led to the removal of 12 misidentified or 
confused sources. The resulting spectral index distribution for 
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Table 5: Euclidean-normalized differential source counts for the GMRT 153 MHz catalogue. 
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Table 6: Binned median spectral indices between the GMRT at 153 MHz and NVSS at 1.4 GHz and WENSS at 327 MHz). 
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a- j 52 is shown in Fig. 10 



Once again there is a bias towards flat- 
ter or inverted spectra"!! low 153 MHz flux densities due to the 
WENSS flux density limit of 18 mJy at 327 MHz. We provide 
upper limits to the spectral indices given the WENSS flux den- 
sity limit for the 576 GMRT sources with that have WENSS flux 
densities below the WENSS detection limit and thus should have 



spectral indices are steeper than the diagonal line in Fig. 10 The 
KM mean spectral index in this case is -0.84 + 0.02 (taking into 
account the upper limits) and is slightly shallower than that ob- 
served between 153 and 1400 MHz. The KM mean spectral in- 
dices measured in 4 flux density bins are also listed in Table [6] 
There is, however, no clear trend with flux density observed, al- 
though there is an indication of a slight flattening of the aver- 
age radio spectrum if the first flux bin is excluded. This may 
be due to the fact that a 3 /E is less robust due to the small fre- 
quency difference and the errors on the individual measurements 
are higher. 

Around 50 per cent of our sources have data at three frequen- 
cies (1400, 327, and 153 MHz), thus we have not attempted to 
fit or locate peaks in the radio spectra. Instead we show a ra- 
dio "colour-colour" plot, Fig. 1 1 comparing the spectral indices 



„1400 



aj^" and cr 153 



3 ^y The line illustrates where the two spectral in- 
dices are equal. Here we have plotted separately bright 153 MHz 
sources, above 0.1 Jy beam -1 , as these sources have smaller er- 
rors on their spectral indices and are not affected by incomplete- 
ness at the other two frequencies (see Figs. |9| and \TU[. In general 



there is a flattening of the average radio spectrum toward lower 
frequencies, as the majority of points fall above the line. It is 
likely that this observed turnover in the spectra at low frequen 
cies is due synchrotron self-absorption. We also plot the distri 



bution of the difference in spectral indices, a 



1400 



„327 
' 153 



, Fig. 
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which shows a mean value of -0.25 for only bright sources anc 
and -0.2 for all sources. 

Finally, we have also compared our source list to VLSS at 
74 MHz ( |Cohen et al.|2007) , again noting that the errors in this 
spectral index will be much greater due to the smaller differ- 
ence in frequency. VLSS has a resolution of 80" so we searched 
for VLSS sources within this radius of each GMRT source. 58 
GMRT sources were matched to 55 VLSS sources. The resulting 
spectral index distribution is shown in Fig. 13 In this case these 
is a bias towards steeper spectra at low 153 MHz flux densities 
due to the VLSS flux density limit of 0.5 Jy at 74 MHz. The KM 
mean spectral index in this case is -0.55 which was calculated 
for sources with GMRT fluxes above 0.5 Jy. 

5. Conclusion 

We have presented the results from a ~ 30 square degree, high 
resolution (25") radio survey at 153 MHz centred on the NOAO 
Bootes field. We have employed the SPAM ionospheric calibra- 
tion scheme to achieve an rms noise in the 7 pointing mosaicked 
image of ~ 2 - 4 mJy beam -1 . The source catalogue contains 
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0.5 - 




153 MHz flux rivl 



Fig. 9: Spectral index, a^™ , distribution of sources matched be- 
tween 1.4 GHz and 153 MHz (grey points). The difference in 
resolution is 45" (NVSS) and 25" (GMRT) and multiple GMRT 
matches to a single NVSS source have been merged into one. 
The vertical line shows 5<T avg , where <T avg is the average rms 
noise in the GMRT mosaic. The diagonal dotted line indicates 
the incompleteness limit due to the sensitivity of NVSS and 
sources with upper limits are plotted as black points along this 
line. The horizontal dashed line shows the KM mean spectral in- 
dex of -0.87 + 0.01 accounting for upper limits. The large black 
triangles show the mean spectral index in 5 logarithmic bins. 
Error bars on individual points are not plotted for clarity, but a 
single bar in the top right indicates the maximum and minumum 
errors in the dataset. 
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Spectral index 

Fig. 11: Comparison between and a 3 ^. The black dashed 
line indicates where the spectral index is the same in both regions 
of the spectrum. Sources with GMRT fluxes above 0.1 Jy beam -1 
are plotted in black and fainter sources are plotted in grey. Error 
bars on individual points are not plotted for clarity, but a single 
error bar in the top left indicates the maximum and minumum 
errors in the dataset. 




153 MHz flux rivl 



Fig. 10: Spectral index, a\y y distribution between 327 MHz and 
153 MHz (grey points). The difference in resolution is 54" x 96" 
(WENSS) and 25" (GMRT) and multiple GMRT matches to 
a single WENSS source have been merged. The vertical line 
shows 5cr avg , where cr avg is the average rms noise in the GMRT 
mosaic. The diagonal dotted line indicates the incompleteness 
limit due to the sensitivity of WENSS. The horizontal dashed 
line shows the KM mean spectral index of -0.84 ± 0.02 which 
takes the upper limits into account. The large black triangles 
show the median spectral index in 4 logarithmic bins. Error bars 
on individual points are not plotted for clarity, but a single bar in 
the top right indicates the maximum and minumum errors in the 
dataset. 



50 



40 




Fig. 12: Comparison between and a™: histogram of 

a i53° ~ a i53 - A£ am ' the histogram for bright GMRT sources is 
plotted in black and for fainter sources in grey. The dashed black 
line shows the mean value of -0.25 and the grey dotted line, the 
mean value of -0.2, for bright and all sources respectively, indi- 
cating that the majority of sources have flattened spectra at low 
frequencies. 



1289 sources between 4.1 mJy and 7.3 Jy detected at 5 times the 
local noise. We estimate the catalogue to be 92 per cent reliable 
and 95 per cent complete to an integrated flux density of 14 mJy. 
The catalogue has been corrected for systematic errors on both 
the astrometry and flux density scales. 

We have analysed the source population by investigating 
the source counts and by identifying counterparts within the 
1.4 GHz NVSS and 327 MHz WENSS surveys and have 
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-1.0 - 



10" 

153 MHz flux flvl 

Fig. 13: Spectral index distribution between 74 MHz and 153 
MHz. The difference in resolution is 80" (VLSS) and the 25" 
(GMRT) and multiple GMRT matches to a single VLSS source 
have been merged into one. The vertical line shows 5cr avg , where 
cr avg is the average rms noise in the GMRT mosaic. The diagonal 
dotted line indicates the incompleteness limit due to the sensi- 
tivity of VLSS. The horizontal dashed line shows the KM mean 
spectral index of -0.55. Error bars on individual points are not 
plotted for clarity, but a single bar in the top right indicates the 
maximum and minumum errors in the dataset. 



computed the spectral index distributions of these sources. 
Understanding the low frequency, low flux source population is 
of particular importance to Epoch of Reionization projects (e.g. 
Ghos h et al.|2012| and references therein) where good models of 
the foregrounds are needed. 

In the near future, this data will be combined with the exist- 
ing multi-wavelength data covering the NOAO Bootes field and 
we will study the properties of radio galaxies as a function of 
various multi-wavelength parameters across a range of cosmic 
time. Further investigation of the spectral indices will be done 
and can be used to identify USS sources as well as high redshift 
gigahertz peaked spectrum (GPS) sources. 
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Fig. A.l: The second brightest source in the 153 MHz catalogue: 
(top) GMRT image and (bottom) FIRST image with GMRT con- 
tours. In both images the GMRT countours are plotted in red 
at intervals of 3cr x [- V3, a/3, VlO, a/30, a/100, . . .] and the 
greyscale goes from 1 <x to 30cr. 



Windhorst, R. A.. Miley, G. K., Owen, F. N., Kron, R. G., & Koo, D. C. 1985. 
ApJ, 289, 494 



Appendix A: Selected Radio Images 



Figures [A2"| shows the 25 brightest sources in the catalogue, ex- 
cluding the second brightest source which is described below 
(see also Fig. A.l I. 



A.1. Note on source J1441 02+3530 

Figure [A~T| shows GMRT postage stamp of the second brightest 
source in the catalogue. Also shown is the FIRST image ( |Becker| 
|et al.|199"5] l of this source which shows that most of the structure 
seen in the GMRT image is in fact real. Only the extension to 
the North-West in the GMRT image has no clear match in the 
FIRST image and may be due to deconvolution errors. 
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Fig. A. 2: The 25 brightest 153 MHz radio sources (exluding J 144102+3530). Countours are plotted in red at intervals of 3<x x 
(- V3, V3, VT6", V30, VTOO, . . .) and the greyscale goes from 1 cr to 30<x. The text in each image lists the local rms noise, the 
source coordinates and total flux, density and the source type ('S' or 'M'). The beamsize is shown in the bottom left corner. 
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Fig. A. 2: Continued. 
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